Quantum correlations in bound-soliton pairs and trains in fiber lasers 
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Quantum correlations in pairs and arrays (trains) of bound solitons modeled by the complex 
Ginzburg-Landau equation (CGLE) are calculated numerically, on the basis of linearized equations 
for quantum fluctuations. We flnd strong correlations between the bound solitons, even though the 
system is dissipative. Some degree of the correlation between the photon-number fluctuations of 
stable bound soliton pairs and trains is attained and saturates after passing a certain distance. The 
saturation of the photon- number correlations is explained by the action of non- conservative terms 
in the CGLE. Photon-number-correlated bound soliton trains offer novel possibilities to produce 
multipartite entangled sources for quantum communication and computation. 
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Solitons in optical fibers are well known for their re- 
markable dynamical properties, at both classical and 
quantum levels. As concerns the latter, the solitons are 
macroscopic optical fields that exhibit quadrature-field 
squeeziiig of quantum fluctuations around the classical 
core as well as amplitude squeezing 0, '^'l , and 

both intra-pulse and inter-pulse correlations |^]. Due to 
the Kerr nonlinearity of the fiber, the quantum fluctua- 
tions about the temporal solitons get squeezed during the 
propagation, i.e., the variance of the perturbed quadra- 
ture field operator around the soliton is smaller than 
in the vacuum state. The nonlinear Schrodinger equa- 
tion (NLSE) is a commonly adopted model for the de- 
scription of classical and quantum dynamics in optical 
fibers. Experimental investigations of quantum proper- 
ties of temporal solitons have shown remarkable agree- 
ment with predictions of the NLSE, provided that loss 
and higher-order effects are negligible iTj, ^ JdJ . 

In the NLSE, adjacent temporal solitons attract or re- 
pel each other, depending on the phase shift between 
them In cases when the potential interaction 

force between the solitons can be balanced by additional 
effects, such as those induced by small loss and gain 
terms in a perturbed cubic [H [l3 or quintic NLSE 
El 111 in El Ei [which actually turns the NLSE equa- 
tion into a complex Ginzburg-Landau equation (CGLE)], 
or the polarization structure of the optical field, described 
by the coupled NLSEs H^,!!!!!!!, bound states of soli- 
tons have been predicted. Recently, formation of stable 
double-, triple-, and multi-soliton bound states [trains, 
in the latter case) has been observed experimentally in 
various passively mode-locked fiber-ring laser systems 
|23, 24, 25], which offers potential applications to optical 
telecommunications. Formation of "soliton crystals" in 
nonlinear fiber rings has been predicted too [261 ]. 



Multiple-pulse generation in the passively mode-lock 
fiber lasers is quite accurately described by the quintic 
CGLE (which is written here in a normalized form). 



iU, + ^Utt + \U\^U = i5U ie\U\^U + il3Uti 
+ i^l\U\'^U -v\U\^U, 



(1) 
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where U is the local amplitude of the electromagnetic 
wave, z is the propagation distance, t is the retarded 
time, and D corresponds to anomalous dispersion (-1-1) 
or normal dispersion (—1). Besides the group- velocity 
dispersion (GVD) and Kerr effect, which are accounted 
for by conservative terms on the left-hand-side of Eq. 
([T|). the equation also includes the quintic correction 
to the Kerr effect, through the coefficient v, and non- 
conservative terms (the coefficients 5, e, and (3 account 
for the linear, cubic, and quintic loss or gain, and spectral 
filtering, respectively). 

In quantum-squeezing experiments, additional noises 
due to the processes other than the GVD and Kerr effect, 
such as the acoustic- wave Brillouin scattering, are un- 
wanted and suppressed, using stable fiber lasers 0. Ac- 
cordingly, the non-conservative terms in the CGLE may 
be superficially considered as detrimental to the observa- 
tion of quantum fiuctuations of fiber solitons. However, 
an accurate analysis of the quantum fiuctuations in the 
CGLE- based model is necessary, and was missing thus 
far, to the best of our knowledge. 

The objective of the present work is to calculate 
quantum fluctuations around bound states of solitons 
in the CGLE by dint of a numerically implemented 
hack-propagation method |27l ]. We find strong quantum- 
perturbation correlations between the bound solitons, de- 
spite the fact the dissipative nature of the model will 
indeed prevent observation of the squeezing of the quan- 
tum fluctuations around the bound solitons in the fiber- 
ring lasers described by the CGLE. Multimode quantum- 
correlation spectra of the bound-soliton pairs show pat- 
terns significantly different from those for two-soliton 
configurations in the conservative NLSE. We also find 
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that a similarity in the photon-number correlations be- 
tween the stable bound-soliton pairs and multi-soliton 
trains. 

Following the known approach to the investigation of 
bound- soliton states |13L Ha . , the corresponding so- 
lution to Eq. Q is sought for in the form 



Uiz,t) 



N 

E 



where Uq is a single-soliton solution, and pj and 9j are 
the coordinates separation and phases of the individual 
solitons. Through the balance between the gain and loss, 
in- phase and out-of-phase bound-soliton pairs may ex- 
ist in the anomalous- GVD regime, which is described by 
Eq. In the case of the normal dispersion, which 

corresponds to the opposite sign in front of Uu on the 
left-hand side of Eq. ^ , strongly chirped solitary pulses 
(which differs them from the classical solitons) and their 
bound states are possible too. 

To evaluate the quantum fluctuations around the 
bound solitons, we replace the classical function U{z, t) in 
Eq. by the quantum-field operator variable, U(z,t), 
which satisfies the equal-coordinate Bosonic commuta- 
tion relations. Next we linearize the equation around 
the classical solution, i.e. U{z,t) = Uq +u{z,t), for a 
state containing a very large number of photons. Then, 
the above-mentioned back-propagation method is used to 
calculate the perturbed quantum fluctuations around the 
bound-soliton states in the full CGLE model. The lin- 
earized equation for the perturbed field operator u(z, t) 



—u{z, t) = -Pi(z, t)u{z, t) + V2{z, t)v){z, t) + n(z, i), (2) 
az 

where Vi and 7^2 are two special operators defined as 
follows, 

-f 3A^|[/o|'^ + 3^^^|C/o|^ 

To satisfy the Bosonic communication relations for the 
perturbed quantum fields u{z,t) and u^{z,t), we also in- 
troduce a zero-mean additional noise operator t) in 
Eq. (jSJ, with satisfy following commutation relations 



[n{z,ti) , h\z',t2)] = 

{-Viiz, h) - Vl{z', t2)}5{z - z')5{ti ~ t2), 
[n{z,ti) , h{z',t2)] = [nHz,ti),h^z',t2)]=0. 

To actually determine the correlation functions of 
n{z,t) and n^{z,t), one has to consider its physical ori- 
gins. In general, h{z,t) = with hi{z,t) 



being the noise operator contributed by the i-th non- 
conservative term in the equation. The commutation re- 
lations of ni{z,t) and are of the same form as 
in the above equations for n{z,t) and n''{z,t), except 
that the differential operators Pi{z,t) and P2{z,t) con- 
tain only the corresponding non-conservative term. For 
simplicity, in our calculation we will introduce the follow- 
ing assumptions: {nl{z,ti)hi{z' ,t2)) — for loss terms 
and {ni{z,ti)hl{z' ,t2)) — for gain terms. With these 
additional assumptions, the correlation functions for each 
t) as well as for the total noise n{z, t) can be calcu- 
lated from their commutation relations. Physically, these 
assumptions are equivalent to assume the reservoirs cor- 
responding to the loss terms are in the ground state and 
the population inversions corresponding to the gain terms 
are in full inversion. The magnitude of the noise level cal- 
culated with these assumptions represents the minimum 
quantum noise that will be introduced with the presence 
of the considered non-conservative terms. For real sys- 
tems, the actual introduced noises will be always larger 
and thus our calculation results here only represent the 
lower bound limit required by the fundamental quantum 
mechanics principles. 

In Fig. ^ we display the comparison of the time- 
domain photon-number correlations for the two- soliton 
configuration in the conservative cubic NLSE model (a), 
and an in- phase two-soliton bound state in the CGLE 
model (b). The correlation coefficients, rjij^ are defined 
through the normally-ordered covariance. 
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An? An? 



(3) 



where An^- is the photon- number fluctuation in the j-th 
slot Atj in the time domain. 



Afij = 



Ati 



dt[Uo{z,t)v){z,t) + UQ{z,t)u{z,t)]. 



Here u{z, t) and Ua{z^ t) are the quantum-field perturba- 
tion and classical unperturbed solution, as defined above, 
and the integral is taken over the given time slot. In the 
NLSE model, there are two isolated patterns for the two- 
solitons configuration, corresponding to the intra-pulse 
correlations of individual solitons, see Fig. Ql^a). How- 
ever, it is obvious from Fig. db) that there is a band 
of strong correlations between the two bound solitons 
in the quantum CGLE model. This strong-correlation 
band can be explained by the interplay between the non- 
linearity, GVD, gain, and loss in the model. The balance 
between these features not only supports the classical 
stable bound state of the soliton, but also causes strong 
correlations between their quantum fluctuations. 

In addition to the time-domain photon-number corre- 
lation pattern, we have also calculated a photon-number 
correlation parameter between the two bound solitons. 
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FIG. 1: A typical pattern of the time-domain photon-number 
correlations, r/ij, for an out-of-phase two- soliton state in the 
cubic NLSE model (a), after 6 normalized propagation dis- 
tance; and for an in-phase two- soliton bound state in the 
CGLE model (b), after 0.5 normalized propagation distance. 



which is defined as 



AN1AN2 :) 
{ANf){ANi) 



Here, AiVi 2 are perturbations of the photon-number op- 
erators of the two sohtons, which are numbered (1,2) ac- 
cording to their position in the time domain. 

Figure |2] shows the evolution of the photon-number 
correlation parameter in the two-soliton bound state. 
Initially, the classical laser statistics (coherent state) is 
assumed for each soliton, without correlation between 
them, C12 « 0. In the course of the evolution, the 
photon-number correlation between the two bound soli- 
tons gradually increases to positive values of C12, and 
eventually it saturates about C12 = 0.36. The inter- 



FIG. 2: The evolution of the photon-number correlation pa- 
rameter, C12, for two bound in-phase. The computations 
were carried out for the following values of parameters in the 
CGLE: D ^1, S = -0.01, e = 1.8, (3 = 0.5, fj. = -0.05, and 
1/ — 0. The inset illustrates the stability of the underlying 
classical solution for the two- soliton bound state, by means 
of contour plots. 



soliton correlation is induced and supported by the in- 
teraction between the solitons. In a conservative system, 
such as the NLSE model, nearly perfect photon-number 
correlations can be established if the interaction distance 
is long enough j28j . In a non-conservative system, such 
as in the CGLE model, the action of the filtering, linear 
and nonlinear gain, and losses lead to the saturation of 
the photon-number correlation parameter. Thus, while 
the large quantum fluctuations in the output bound- soli- 
ton pair will eclipse any squeezing effect, the correlated 
fluctuations between the bound soliton are predicted to 
be observable. 

The approach elaborated here for the study of the 
quantum-noise correlations between two bound solitons 
can be easily extended to multi-soliton bound states (soli- 
ton trains). Figure shows the photon- number correla- 
tion parameters, Cy , in a train of four equal-separated 
in-phase bound solitons. Again, the photon-number fluc- 
tuations are initially uncorrelated between the solitons. 
As could be expected, we find that soliton pairs with 
equal separations have practically identical correlation 
coefficients, i.e. , 6*12(2:) « 6*23(2;) w (734(2;), and 
Ci3(2) ~ (724(2). Obviously, when the interaction be- 
tween the soliton trains is stronger, as the separation 
between them is smaller, the values of the correlation pa- 
rameter are larger. Note that the correlation coefficient 
for the most separated pair in the train, C14, grows very 
slowly with z. Eventually, all of the curves of Cij{z) sat- 
urate due to the dissipative effects in the CGLE model. 

In conclusion, in this paper we have extended the 
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FIG. 3: The evolution of the photon-number correlation pa- 
rameters, Cij, for a bound complex (train) of four in-phase 
solitons in the CGLE model. The soliton's number (1, 2, 3, 4) 
runs from left to right. Parameters are the same as in the 
case shown in Fig. The inset demonstrates the stability 
of the underlying classical solution for the four-soliton train. 



concept of quantum fluctuation around classical opti- 
cal solitons to the non-conservative model based on the 
CGLE (complex Ginzburg-Landau equation) with the 
cubic-quintic nonlinearity. Applying the known back- 
propagation method to the linearized equations that gov- 
ern the evolution of the fluctuations, we have numerically 
calculated the photon-number correlations and the effec- 
tive correlation coefficient for pairs of bound solitons, as 
well as for multi-soliton bound complexes (trains). We 
have demonstrated that, unlike the two-soliton configu- 
ration in the conservative NLSE model, there is a band 
of strong quantum correlation in the bound-soliton pair. 
While the dissipative effects in the CGLE model will 
totally suppress the generation of squeezed states from 
bound solitons, as one might expect, there still exists 
a certain degree of correlations between photon-number 
fluctuations around the stable bound-soliton pairs and 
trains. Recently, experimental progress in the study of 
various quantum propert ies of solitons in optical flbers 
has been reported [l^, yfl IM, |32| , which opens the way to 
observe effects predicted here. Besides that, the photon- 
number-correlated soliton pairs and trains, predicted in 
this work, may offer new possibilities to generate mul- 
tipartite entangled sources for applications to quantum 
communications and computation. 
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